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ABSTRACT 

In  this  paper,  a  complete  description  of  dynamic 
modeling  of  proton  exchange  membrane  fuel  cell 
(PEMFC)  and  solid  oxide  fuel  cell  (SOFC)  is  given. 
For  modeling  of  fuel  cell  for  software  based  study, 
calculation  of  all  voltage  drops  during  chemical 
reaction  of  fuel  is  required.  Additionally  a  flow  chart 
of  fuel  cell  output  voltage  calculation  is  also 
explained  which  includes  fuel  cell  voltage,  double 
layer  charging  effect,  thermodynamic  response  and 
terminal  voltage  of  fuel  cell.  By  using  ac  to  dc 
converter,  the  fuel  cell  power  can  be  connected  to 
load  or  grid.  Based  on  this  study,  a  mathematical 
model  of  fuel  cell  is  developed  in  simulink  software 
MATLAB  to  obtain  output  characteristic  of  fuel  cell. 
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I.  INTRODUCTION 

Fuel  cells  are  static  energy  conversion  devices  which 
convert  the  chemical  energy  of  fuel  directly  into  dc 
electrical  energy.  The  physical  structure  of  fuel  cell 
consists  of  two  porous  electrodes  and  an  electrolyte 
layer  in  middle. 

Fuel  cells  are  good  energy  sources  to  provide  reliable 
power  at  steady  state  but  they  cannot  respond  to 
electrical  load  transients  as  fast  as  desired.  This  is 
mainly  due  to  their  slow  internal  electrochemical  and 
thermodynamic  responses  [1].  Fuel  cell  models  are 
needed  to  predict  and  evaluate  their  steady-state  and 
dynamic  responses  in  various  applications,  such  as 
fuel  cell  vehicles  and  stand-alone  and  grid-connected 
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fuel  cell  power  generation  systems.  Dynamic  fuel  cell 
models  are  also  needed  for  predicting  and  evaluating 
fuel  cell  performance  under  different  transient 
conditions  such  as  electrical  faults  at  the  fuel  cell 
terminals.  In  these  conditions,  output  voltage  and 
power  of  fuel  cell  are  affected  by  physical  properties 
of  fuel  cell  such  as  voltage  drop  inside  fuel  cell, 
capacitance  of  double  layer  charge  effect,  material 
conversion  etc.  Dynamic  models  of  fuel  cells  are  also 
needed  for  designing  controllers  to  control  fuel  cell 
terminal  voltage  as  desired.  In  the  applications  where 
fuel  cells  are  connected  to  grid  or  operated  in 
standalone  system,  it  is  possible  to  control  fuel  cell 
performance  by  controlling  power  electronic  devices. 

The  fuel  cell  output  voltage  is  less  than  internal 
voltage  developed  inside  fuel  cell  due  to  several 
voltage  drops.  The  ohmic  voltage  drops  associated 
with  electrodes  and  electrolyte  at  certain  operating 
point  can  be  considered  as  linear  function  of  thickness 
of  those  segments.  The  activation  voltage  drop  is  due 
to  energy  barrier.  The  concentration  voltage  drop  is 
due  to  current  density  different. 

There  is  growing  interest  to  use  hydrogen  and  fuel 
cells  to  power  stationary  and  mobile  (transportation) 
applications  due  to  their  high  efficiency,  modularity, 
environmental  friendliness  and  compact  design  [2], 
An  individual  fuel  cell  produces  less  than  a  volt  of 
electric  potential.  A  large  number  of  cells  are  stacked 
on  top  of  each  other  and  connected  in  series  (with 
bipolar  connects)  to  produce  higher  voltages  [3]. 
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II.  Dynamic  modeling  of  PEMFC 

The  output  voltage  of  PEMFC  is  given  by  Nemst 
equation 

RT 

Ecell  =  Eo.cell  +  —In  [PniiPoz)03] 

Where, 

Ecell  =  Potential  of  each  cell 
Eo,cell  =  Reference  potential 
F  =  Faraday’s  constant 
R  =  Gas  constant 
T  =  Fuel  cell  temperature 
Ph2  =  Pressure  of  hydrogen 
P02  =  Pressure  of  oxygen. 

The  overall  effect  of  the  fuel  and  oxidant  delay  is 
represented  by  a  voltage 

t 

Ed, cell -=  —  i(t)aer] 

Where, 

D  is  convolution  operator  and  et/T  is  time  lag  function. 
So  total  voltage  generated  is  given  by, 

RT 

Ecell  ~  Eo.cell  T  ^p  ^[^2(^02)°  5]  —  Ed  ceu 

Under  normal  operating  condition,  the  fuel  cell  output 
voltage  is  less  than  voltage  generated  due  to  ohmic 
voltage  drop,  activation  voltage  drop  and 
concentration  voltage  drop. 

Eceii  —  Eceu  EaC(  Vohm  ~  Krone 

And  Vout  =  Ncell  x  Vcell 

Where  Ncell  is  number  of  fuel  cell  stack 

A.  Activation  voltage  drop: 

The  activation  voltage  drop  is  a  function  of  the 

PEMFC  current  and  temperature.  Vact  can  be 

described  as  sum  of  Vacti  and  Vact2[l]- 

Vact  =  Vacti  +  Vact2 

Vact  =  [r|o  +  (T-298)*a]  +  [T*b*ln(i)] 

Whereto, a  and  b  are  empirical  constants. 

Vacti  is  the  voltage  affected  only  by  the  fuel  cell 
internal  temperature  and  Vact2  is  both  current  and 
temperature  dependent.The  equivalent  resistance  of 
activation  voltage  drop  is  ratio  of  Vact2  and  fuel  cell 
current. 

Pact  —  Eact2^E 

B.  Ohmic  voltage  drop: 

The  ohmic  resistance  of  PEMFC  consists  of  resistance 
of  polymer  membrane  and  resistance  of  electrodes. 

V ohm  —  V memb  V ele 


C.  Concentration  voltage  drop: 

This  voltage  drop  depends  on  current  density.  At 
higher  current  density,  the  reaction  of  chemicals 
becomes  slow  and  it  is  due  to  concentration  voltage 
drop.  The  equation  of  this  voltage  drop  is  given  as, 


Where,  z  =  numbers  of  electrons 
hm= fuel  cell  limit  current 

D.  Double  layer  charge  effect: 

As  the  electrons  flow  from  cathode  to  anode  through 
external  load,  the  electrodes  act  ac  two  layers  of 
capacitor  having  the  charge  and  the  membrane 
between  both  acts  as  dielectric  between  both 
electrodes  forming  a  capacitor. 

/  dVc\ 

Vc  =  ^1  —  C  x  ( Ract  +  Rconc ) 

By  considering  all  these  voltage  drops  and  double 
layer  charge  effect,  the  output  voltage  is  given  by, 
Vout  =  E  —  Vact  —  Vc  —  Vohm 

III.  Dynamic  modeling  of  SOFC 

An  SOFC  based  Distributed  generation  System  is 
normally  interfaced  with  the  utility  network  through  a 
set  of  power  electronics  devices.  The  interface  is  very 
important  as  it  affects  the  operation  of  the  fuel  cell 
system  and  the  power  grid.  Various  control  schemes 
haven  proposed  in  the  recent  work  to  interface 
different  energy  sources  to  utility  grid.  Nemst 
equation  to  give  output  voltage  of  Solid  Oxide  fuel 
cell  is  given  by, 

RT  (Pin)2  P02 

Ecell  =  E0  ceu  +  —  In  [y  )  °  ] 

4F  (A) 

Where  Eo,Ceii  is  temperature  dependent  term. 

Similar  to  PEMFC  output  voltage,  same  equation  of 
Vceii  is  used  to  calculate  actual  voltage  generated  in 
fuel  cell  i.e.  it  also  includes  activation,  concentration 
and  ohmic  voltage  drops. 


A.  Activation  voltage  drop: 

This  voltage  drop  is  due  to  activation  of  energy 
barrier  when  there  is  a  chemical  reaction  inside  a  fuel 
cell  at  lower  currents.  The  activation  voltage  drop  for 
SOFC  can  be  estimated  by  following  equation. 

RT  i 


Vact  = 


zfF  ln  W 


This  equation  is  called  Tafel  equation  in  which  the 
value  of  voltage  drop  approaches  infinite  if  current  is 
zero  which  is  undesirable.  The  significance  of  R,  T,  z 
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and  F  is  same  as  in  PEMFC.  The  value  of  p  is  0.5  for 
solid  oxidefuel  cell  [5]. 

B.  Ohmic  voltage  drop: 

The  resistances  due  to  electrodes,  electrolyte  and 
interconnection  of  several  fuel  cells  are  responsible 
for  total  Ohmic  voltage  drop 

Vohm  =  Velectrode  +  Velectrolyte 
+  Vinter  connection 


C.  Concentration  voltage  drop: 

For  higher  current  densities  and  due  to  slow  diffusion 
process,  the  reactants  cannot  be  transported  to 
reaction  site  as  fast  as  needed  and  this  causes  a 
significant  voltage  drop  called  Concentration  voltage 
drop.  This  voltage  drop  can  be  obtained  by  calculation 
voltage  drops  at  cathode  and  anode  terminals. 

RT 

Vconc  =  -—In 
4  F 


XPH2)2Po2 

RT 

- In 

XPCH2)2Po2 

(PaH2oY  \ 

4  F 

(PCH2oY  \ 

D.  Double  layer  charge  effect: 

Similar  to  PEMFC,  a  double  layer  charge  effect  is 
present  in  SOFC  in  same  manner  and  can  be 
estimated  by  same  equation  as  in  PEMFC.  The 
difference  is  that,  in  SOFC,  the  electrodes  are  porous 
and  it  makes  the  value  of  capacitance  to  be  higher 
comparatively. 


After  considering  the  voltage  drops  and  capacitor 
effect,  the  output  voltage  of  SOFC  can  be  calculated 
as  below. 

Vout  =  Ecell  —  Vact  —  Vconc  —  Vohm  —  Vc 


IV.  Flow  chart  of  fuel  cells  output  voltage 
calculation 

For  the  mathematical  modeling  of  fuel  cell,  a  step  by 
step  approach  in  a  form  of  flow  chart  is  shown  below. 
This  includes  input  parameters  and  calculation  of 
different  voltage  drops. 

As  shown  in  figure,  the  input  parameters  of  fuel  cell 
are  Pressure  at  anode  and  cathode  terminal  (Pa  and 
Pc),  room  temperature  (Tr00m),  Initial  temperature 
(Tinitiai),  Air  temperature  (Tair),  fuel  temperature  (T  iuei), 
hydrogen  flow  rate  (MH2),  water  vapour  flow 
rate(MH20)-  These  parameters  are  applied  to  Nemst 
equation  to  calculate  the  fuel  cell  internally  generated 
voltage.  Here  the  material  conservation  and 
thermodynamic  effect  of  fuel  cell  is  neglected.  So  the 
voltage  drops  can  be  calculated  as  per  the  equations 
and  voltage  drop  due  to  capacitor  effect  in  fuel  cell 
can  be  calculated.  The  addition  of  all  calculated 
voltage  drops  when  subtracted  from  voltage  estimated 


by  Nemst  equation,  gives  terminal  voltage  of  fuel  cell. 
The  equations  and  parameters  selected  for  modelling 
are  important  for  analysis  and  control  of  fuel  cell  [6], 


Fig  1:  Flow  chart  of  fuel  cell  output  voltage 
Calculation 

V.  Simulink  model  of  fuel  cells 

For  dynamic  modelling  of  fuel  cells,  mathematical 
equations  of  voltage  drops,  Nemst  equation  are 
modelled  in  Simulink  software  MATLAB  to  get  the 
output  voltage.  The  mathematical  model  of  PEMFC  is 
shown  in  Fig  2.  The  standard  values  of  parameters  for 
calculation  of  Fuel  cell  output  voltage  is  given  below. 
In  this  model,  the  thermodynamic  effect  of  fuel  cell  is 
neglected. 

Table  I.  PEMFC  parameters 


Parameter 

Value 

Gas  Constant  R 

8.3143  J/mol*K 

Faraday’s  Constant  F 

96485  J 

Pressure  of  anode 

1  Atm 

Pressure  of  cathode 

1  Atm 

Absolute  Temperature 

298.16  K 

Reference  Temperature 

273.15  K 

Ilim 

25.1  Amp 

Number  of  cells 

24 

Double  layer  capacitance 

1  F 

Heat  transfer  coefficient 

10  W/K 
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Fig  2.  Mathematical  model  of  PEMFC 


Considering  above  parameters,  the  model  of  PEMFC 
is  developed  and  voltage-current  characteristic  is 
shown  in  fig  3. 


I.(r 


-  1 


is - 

W - 1 - ^ - 1 - 1 - 1 - 1 - 1 - 1 

0  1  2  3  4  5  6  7  8 

current(A) 

Fig  3.  PEMFC  V-I  Characteristic 


In  V-I  characteristic,  we  can  see  that  the  output 
voltage  is  going  down  because  of  three  voltage  drops. 
Initially  at  lower  current  values  activation  and 
concentration  voltage  drop  is  higher  and  than  ohmic 
drop  causes  further  reduction  in  output  voltage. 


A  mathematical  model  of  SOFC  can  also  be 
implemented  by  simulating  Nernst  equation  and 
voltage  drops  to  obtain  output  voltage.  The 
parameters  considered  for  calculation  are  given 
below. 

Table  II.  SOFC  parameters 


Parameter 

Value 

Gas  Constant  R 

8.3143  J/mol*K 

Faraday’s  Constant  F 

96485  J 

Ideal  standard  cell 

potential 

1.18  V 

Foad  Resistance 

10  Q 

Absolute  Temperature 

1273  K 

Reference  Temperature 

273.15  K 

linitial 

100  Amp 

Molar  Constant  H2 

8.43  e-4  K  mol/(s  atm) 

Molar  Constant  O2 

2.52  e-4  K  mol/(s  atm) 

Molar  Constant  H20 

2.81  e-4  K  mol/(s  atm) 

For  this  model  of  SOFC,  with  above  parameters,  the 
output  voltage  of  1.26V  is  obtained.  The  advantage 
with  SOFC  is  that,  as  a  working  temperature  is  high 
(600°C  -  1000°C)  it  can  be  used  with  combined  heat 
and  temperature  so  as  to  give  higher  efficiency. 
According  to  the  load  profile,  the  feasible  fuel  cell 
capacity  can  be  determined,  whereas  a  diesel 
generator  should  be  operating  at  the  rated  power  as 
much  as  possible  [4], 
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VI.  CONCLUSION 

A  software  based  dynamic  modeling  of  fuel  cell  can 
be  used  to  predict  its  performance  for  different 
operating  conditions  and  parameters,  by  the  V-I 
characteristic.  The  voltage  drop  inside  the  fuel  can 
also  be  estimated  by  this  model  and  initiative  can  be 
taken  to  reduce  the  voltage  drop  to  increase  the  output 
voltage  and  efficiency.  The  dynamic  behavior  of 
PEMFC  and  SOFC  can  also  be  assessed  by  changing 
the  pressure,  temperature  or  flow  rates  of  fuel  and 
oxidant. 
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